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Abstract 

The dynamics of a membrane is a coupled system comprising a moving elastic surface 
and an incompressible membrane fluid. We will consider a reduced elastic surface model, 
which involves the evolution equations of the moving surface, the dynamic equations of 
the two-dimensional fluid, and the incompressible equation, all of which operate within 
a curved geometry. In this paper, we prove the local existence and uniqueness of the 
solution to the reduced elastic surface model by reformulating the model into a new 
system in the isothermal coordinates. One major difficulty is that of constructing an 
appropriate iterative scheme such that the limit system is consistent with the original 
system. 

1 Introduction 

This paper is concerned with the hydrodynamics on the moving surface of bio-membrane, 
which as the outerwear of living cehs and organelles plays an important role in the life 
process. Consisting of lipids, proteins and carbohydrates, the structures and properties of 
bio-membrane are very complex. In general, bio-membrane can be viewed as a 2-dimensional 
fluid surface consisting of a lipid bilayer, as the lipid molecules can move freely on the surface 
but cannot escape from it. The fluid is viscous and can be viewed as incompressible because 
it typically has a large tensile module. Moreover, this 2-dimensional fluid surface is bend- 
resistent. Hence, it tends to minimize the Helfrich energy under the fixed area condition 
(guaranteed by the incompressible condition) 



Eh = J {ci{H - Bf - C2K)da, (1.1) 



where H and K are the mean curvature and the Gaussian curvature, respectively, B is the 
spontaneous curvature that reflects the initial or intrinsic curvature of the membrane, ci and 
C2 are the elastic coefficients, and dcr is the area form of the surface [6]. When C2 is uniform 
on the membrane, J Kda is a constant determined by the topology of the membrane. When 
S = 0, Eh is called the Willmore energy in geometry. A number of studies based on Helfrich's 
bending energy model explore the mechanics of bio-membrane, for example, see [221 [3| [H]. 
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During the past several decades, membrane dynamics have received considerable atten- 
tion. Researchers from different fields have developed several models with/without the sur- 
rounding fluid to study the behaviors of the membrane. For the models without surrounding 
fluid, see [251 |26l Ell [Ml E], and for the models with surrounding fluid, see [171 ITSl IT3l [T6] . 

Waxman [25] may have been the first to study the dynamics of bend-resistant bio- 
membrane using a model without surrounding fluid and in which the incompressibility, 
bend-resistance, and viscosity effects are all considered. However, Waxman's model does 
not preserve the energy dissipation law. In jlOj . Hu-Zhang-E introduced a director field to 
represent the direction of lipid molecules at every material point and developed an elastic 
energy model based on the Frank energy of the smectic liquid crystal. When the director is 
constrained to the normal of the surface, they obtain a reduced elastic surface model, that 
is very close to Waxman's model, but adds one term to the in-plane stresses whereby the 
model satisfies a natural energy dissipation law. In the elastic surface model, the dynamics of 
the membrane involves the evolution equations of the moving surface, the dynamic equations 
of the two-dimensional fluid, and the incompressible equation, all of which operate within a 
curved geometry. 

For a surface membrane F = R(ui,M2,t), we denote by the tangent vectors of F, n 
the unit normal vector, (aa/3)i<o,/3<2 the covariant metric tensor, Ar the Lapalace-Beltrami 
operator, K the Gaussian curvature, and H the mean curvature. In the simple case, the 
reduced elastic surface model takes the following form: 

^ = v{ui,U2,t), 

^ = (-na"/^a,),/3 + 2eo(5"''a,),^ - l(ArF + 2H{H^ - K))n, (1.2) 
Vr • V = 0. 

Here, v is the velocity of the fluid, H is the surface pressure, S'"'^ is the rate of the surface 
strain, and the constant eo > is the shear viscosity. The notation ()^^ denotes the covariant 
derivative. The first term on the right-hand side of the second equation is induced by the 
incompressible condition Vp-v = 0, and the surface pressure H can be viewed as a Lagrangian 
multiplier; the second term describes the viscosity of the fluid on the surface; the third term 
is the elastic stress induced by the Helfrich bending energy (jl.ip with B = 0. Please see 
Section 2 or [10] for more detail. 

When the interaction with bulk fluid is considered, Hu-Zhang-E [10] also derived the 
incompressible membrane-fluid coupling system in the form 

U( + u • Vu = — Vp + z^Au, in fi, 

V • u = 0, in Q, 

[-pi + r] • n = F, on F, (1.3) 

[u] = 0, on F, 

Vr • u = 0, on F, 

where r = i^(Vu + Vu^) is the stress of the bulk fluid, F is given by the right-hand side 
of the second equation of (jl.2p . Q is the fluid domain, F is the time-dependent surface of 
the membrane included in Q, and [•] denotes the jump across the membrane. In a recent 
review paper [16], a similar model was derived via the direct variational method. Compared 
with the classical free boundary problem of the Navier-Stokes equations, the main difference 
is that the system (II. 3p contains two unknown pressures: the pressure p of the surrounding 
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fluid and the pressure IT of tlie membrane defined on the surface, where 11 is determined by 
the incompressible condition Vr • u = 0. Due to the coupling between p and H, solving the 
membrane-fluid coupling system ()1.3p is still challenge, both mathematically and numerically. 
In some specific case (e.g., when the velocity of the surrounding fluid is small), the main 
influence of the bulk fluid is to maintain the enclosed volume of the membrane. In such cases 
for simplicity, it can be replaced by introducing osmotic pressure. Moreover, although the 
reduced model (jl.2p neglects the fluid interaction, numerical simulation [8] also convinces 
us that this model can be used to reconstruct some important physical processes, such as 
exocytosis and endocytosis. 

To our knowledge, few mathematical results such as the well-posedness for the fluid bio- 
membrane dynamics are available. In [3], Cheng-Coutand-Shkoller studied the bulk fluid 
interacting with a membrane considered a nonlinear elastic bio-fluid shell and modeled by the 
nonlinear Saint Venant-Kirchhoff constitutive law, where the membrane is compressible and 
the surface fluid is inviscid. In |9j, Hu-Song-Zhang proved the local existence and uniqueness 
of ()1.2p for a simplified case when the membrane is cylindrical. In this case, the membrane is 
similar to a 1-D incompressible string such that the fluid vanishes. With the introduction of 
the arc length parameter and the tangent angle of the curve, the system is transformed into 
a fourth-order wave equation for the tangent angle a coupled with an elliptic equation: 

f att =91 + ^TsOs + Toss -{B + as)sss + olI{B + as)s, 
\ -Tss + a^T = 52 + a? + 2{B + a,)ssas + {B + Qs)sass, 

where gi,g2, and B are the given smooth functions. 

The purpose of this paper is to prove the local well-posedness of the system ()1.2p . This is 
also a key step toward understanding and solving the membrane-fluid coupling system (jl.Sp . 
Our result is stated as follows. 

Theorem 1.1 Let s = 2k for some integer k > 3. Assume that the initial velocity vq G H^~^ 
and the initial closed surface Rq S H^~^^. There exists T > such that the system has 
a unique solution (v(t),R(i)) on [0, T] satisfying 

V G C{[0,T];H'-^), R G C{[0,T]- H'+^). 

Remark 1.2 The regularity we imposed on the initial data should not be optimal. To simplify 
the analysis, we will work in a functional space with high regularity. 

System ()1.2p is a coupled system of parabolic, hyperbolic, and elliptic equations. The 
evolution equations of the tangential velocities are parabolic, the evolution equations of the 
normal velocity and the mean curvature constitute a hyperbolic system, and the pressure 
satisfies an elliptic equation, see (|3.23p - (l3.28j) . Because the surface is moving, it seems natural 
to solve (II. 2p in the framework of Lagrangian coordinates. However, some essential difficulties 
will arise. Let us explain it in what follows. 

Assume that the initial velocity vq G H^~^ and the initial surface Rq G H^~^^. Because 
the tangential velocity satisfies the parabolic equation, and the normal velocity f " and the 
mean curvature H together satisfy the hyperbolic system, it seems natural to expect v" to 
belong to L^(0, T; i/**), and {v"',H) to belong to L°°{0,T; H^~^). However, these estimates 
depend on the regularity of the metric of the surface. Hence, we have to recover the 
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regularity of the metric from (v°' ,v'^ , H) in order to close the energy estimates. In the 
Lagrangian coordinates, we have 

Rj = v(ni,n2,t), 

which tells us that R € L°^(0, T; H^~^) by the estimate for the velocity. Hence the metric has 
only H^~^ regularity (a loss of two derivatives). Maybe, one wants to use the regularity of the 
mean curvature to gain the regularity of R(Note that formally, H^~^ regularity of the mean 
curvature suggests that the free surface has H^^^ regularity). However, we cannot expect 
R to have more regularity in the Lagrangian coordinates, see the example and argument of 
Section 5 in [20]. 

Another way to solve the system is to represent the moving surface locally by 2:3 = 
g{xi,X2,t), where g satisfies the following hyperbolic equation 



qtt + Arfdivrf — =)] 



lower-order terms. 



However, if we make an energy estimate for this equation, the estimate is also not closed, 
since the lower-order terms contain the third-order derivative of g, which cannot be controlled 
by the main part. 

Motivated by [T], we will use the isothermal coordinates to re-parameterize the surface. 
There are two main advantages adopting the isothermal coordinate: (1) we can gain two more 
regularities for the surface from the regularity of the mean curvature, and (2) the coefficients 
of the first fundamental form have the same regularity as the surface. Indeed, there are the 
following important relations between the surface T = R(ni,M2), the first fundamental form 
E, and the mean curvature H when (ui, U2) is taken as the isothermal coordinate of F: 

AR = 2EHii, AE = 2{du,du,B.-d^,d^,R-dl^R-dl^R). (1.4) 

Here A = d^_^ + d'^^, and n is the unit normal of F. 

In general, it is difficult to construct an approximate system preserving the isothermal 
relation. As the solution of the approximate system does not satisfy the important geometric 
relation (jl.4p . there will also be derivative loss once we make the energy estimates for the 
approximate system. To overcome this difficulty, we incorporate the relation (jl.4p into our 
iterative scheme. However, this produces another very troubling problem — one that arises 
for the most part from the construction of the iterative scheme and relates to the equivalence 
of the two systems. The problem is this: we do not know and need to establish whether the 
limit system is equivalent to the original system, and proving such equivalence involves very 
complicated geometric calculations. 

This paper is organized as follows. In the next section, we review some formulae for the 
evolving surfaces and introduce the reduced elastic surface model. In Section 3, we derive an 
equivalent system in the isothermal coordinate by decomposing the velocity into tangential 
and normal components. Section 4 is devoted to studying the linearized system. In Section 
5, we prove our main results, including the construction of the iteration scheme, nonlinear 
estimates, the convergence of the iteration procedure, and the consistency between the limit 
system and the original system. 
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2 The elastic model of an incompressible fluid membrane 



In this section, we provide a short derivation of the dynamic model of an incompressible 
elastic fluid membrane in three-dimensional space. We refer to [10) for more details. 



2.1 Geometric tensors and their evolution equations 

For a surface membrane F = R({i, t) with a curve coordinate u = {u^^u^), we can get the 
Frenet coordinate system of the surface. Namely, the tangent vectors and the unit normal 
vector n are given by 

d'R ai X a2 

ou'^ |ai X a2| 

The covariant metric tensor {aaj3)i<a,p<2 is defined as 

ao/3 = aQ, • a/3. 

We denote its inverse by {a'^^)i<.a,p<2-, which can be used to raise or lower the indices of the 
vectors and tensors. For example. 

The surface Christoffel symbols F^^ = F^^ and the curvature tensor 6^^ = are given by 
the Gauss- Weingarten-Codazzi equation: 



Here we use a comma followed by a lowercase Greek subscript to denote the covariant deriva- 
tives based on the metric tensor a^^, that is, 

Qt.,. = -i^ + T.n.Qt. - Y.n.Qt.- (2-1) 

For example, we have 

Oa/3,7 - -g-f - a-yOSfi " i f}^bo,5- 

Thus we can rewrite the Gauss- Weigarten-Codarzzi equation as 

aa,/3 = baf^n, Ha = —ba^p, bal3,'y = ba^^p. (2.2) 

The mean curvature H and the Gaussian curvature K of the surface are given by 



H = \bl, K = \{AH'-bpi). 



2 2 

In the following, let us derive the evolution equations of the geometric tensors. For this 
purpose, we denote by v(n, t) the velocity of the surface given by 

^{u,t) = —^, (2.3) 
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and we decompose it into 



V = v"aa + v'^n. 



Using (j2.2p . it is easy to find that 

dsia d dH d dH dv 



dt dtdu°' du°' dt du°' 

= {vl-v^bl)si, + {v% + vn^^)n. (2.4) 



As n • aa = 0, we get by (f2^ that 

dn dsia 



a, 



dt " dt 



which together with the fact of ^ • n = imphes that 



dt 

^ = -(^/36^ + a-V^)a,. (2.5) 
The evolution equation of the metric tensor is given by 

—Qf = -Qf ■^f} + -Qf^a = {Va,f^ + Vi3,a) - 2v b^f^. (2.6) 

Differentiating the identity a'^^apr^ = 5'^ with respect to t, we get by (j2.6p that 

= -a-^a^^y^^ + v,^^) + 2^;«6-/^. (2.7) 
And differentiating bap = — a^, • n^^, we get by (|2.2p - (j2.5p that 

%^ = {v:ap - v'^blb.p) + {v^b^^ + vlb.p) + vn^p^,. (2.8) 



Due to 2H = a°'^baf3, we get by ([22]) and ([ZH]) that 

2^ = a'^^v^^ + v^b-pbi + 2v-H,^. (2.9) 
2.2 The derivation of the elastic surface model 

In this subsection, we choose u = {v},v?) as the Lagrangian coordinate of the moving fluid 
surface. In this coordinate system, the velocity of the fluid on the surface is equal to the 
velocity of the fluid v given by (j2.3p . The Helfrich bending elastic energy (6l [7] is 

Eh = Cf^\B^p - b^p){B^s - ^5)d5, (2.10) 
where B^^p is the spontaneous curvature tensor, and the fourth-order tensor C'^^^^ is given 

by 
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where ki and ei are positive elastic coefficients and ki > ei. 

As the membrane is a two-dimensional incompressible fluid, we have 

Vr • V = 0, 

which is equivalent to 

v% = 2Hv'^. (2.11) 

Then by applying the principle of virtual work, we obtain elastic stresses. For isotropic 
Newtonian membrane fluids, the dynamical equation of the membrane is 

^|^=f+(r"%),„ + (g-n),,, (2.12) 

where g is the membrane fluid density. The in-plane stress tensor T"-^ and transverse shear 
stress g° are given by 

= Mf, 

where IT is the surface pressure (tension), and the rate of the surface strain is given by 

def 1 daai3 1 . , \ n. 

From (j2.1ip . it is easy to see that = 0, and thus J"^ = leoS'^^. Furthermore, the above 
equations have the following energy dissipation relation: 

^^^(^Eh + I^Q\^r\''ds') =- j^C^^'<^S^pS^sdS = -2eo j^S^^S^pdS. (2.13) 

If the function B{u) = B[u)aaj3 with B independent of the time then ()2.10p can be reduced 
to 



Eh = j iki{H - Bf + 4:ei{H^ - K)dS, 



and the velocity equation (j2.12p can be reduced to the following form 

= {-Ua''^Sic.),p + 2eo{S''^Sia),p-^kiHa''^B^paia 

+2ki{ArB - 2BK)n - 2{ki + ei)(^ArH + 2H{H^ - K)^n. (2.14) 

Here Ar denotes the Lapalace-Beltrami operator on the surface F, and K is the Gaussian 
curvature. We refer to the appendix for the derivations of (I2.13P and (|2.14p . 

In this paper, we only consider the simple case with B{u) = 0. By the rescaling argument, 
we can set g = 1, 4(A;i -|- ei) = 1. Thus, we obtain 

f f = (-na^/^a,),;? + 2eo(S"^a,),;3 - i (ArF + 2H{H^ - K)) n, 
I v'i = 2Hv''. 
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3 New formulation of the system 



Motivated by [T], we will reformulate (j2.15p into a new system in the isothermal coordinates. 
That is, we choose a coordinate x = {x^,x'^) such that the metric tensor satisfies 



an = 022, ai2 = 021 = 0. 



(3.1) 



As the tangential velocity of the surface only serves to reparameterize the surface, in the 
following we choose them such that 



dan da22 dai2 



dt 



dt ' dt 



0. 



(3.2) 



As a result, if it holds for the initial surface the relation (jS.ip will be preserved for any time 
t. 



3.1 Elliptic system for the tangential velocity of the surface 



In the sequel, for convenience, we denote fa = '§^ifi3 = ) a-iid whereas (•)^q, (or (•),/?) 



denotes the covariant derivative with respect to (or x^). The unit tangent vector and the 
unit normal vector of the surface are given respectively by 



Ro 

iR-n 



n 



'/3 



Rn. X Rfi 



(3.3) 



We denote E, F, and G by the coefficients of the first fundamental form, and L, M, and N 
by the coefficients of the second fundamental form. Namely, 



E 
L 



Ra • ; F — • R/3 , G — R^ • R/3 , 



Rn 



n, 



M = R, 



■Q/3 • 



n. 



iV = R 



■/3/3 • 



n. 



In the isothermal coordinates, we have E 
calculated as follows: 



11 

12 



12 
J- 21 



J" 21 

p2 

22 



G,F = 0. The Christoffel symbols can be 

(3.4) 



J" 22 
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Eg 

2E' 
2S- 





t2 = 








tl = 








n = 




riQ, 




n = 


-t2. 


riQ, 




n = 




n/3 




n = 


-t2. 


n/3 



And the following identities can be verified easily: 



R;3 
IR-/3I 

|Ra I 
Rg 
|Rg| 

Ra 
IR/3I ■ 

Ra 
|Rg| 

|R/3l ■ 



riQ, 
n/3 



Rq 


Rg^Rg 


Ep 


|Ra| 


E ~ 


2E-> 


R/3 


Rg/3-R^ 


Ed 




E ~ 


2E^ 














M 






Ve' 






M 






Ve' 






N 






~ Ve- 







(3.5) 



For a given normal velocity U^{x^ , x'^ ,t), we assume that the evolution of the surface is 
determined by 

^^(^^'Ku-n + W,t' + W2t''^'^ 



dt 



w. 



(3.6) 
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Then it follows from (j3.5p that 



def 



def 

Consequently, we obtain 



Ain + ^oit^ +^3t^ (3.7) 
Aan + ^4t^ + Ao2t^ (3.8) 



Et = 2Rc,t • Ra = 2\/^Rc,t • = 2VEA01, 
Gt = 2Tipt ■ R/3 = 2VEKpt ■ = 2VEA02, 
Ft = VE{K^t • t2 + R^i . ti) = Ve{A3 + ^4) 
Now the relation (|3.2p is equivalent to 

{E - G)t = 0, Ft = 0, 

which implies that 

f wA _(W2\ _ t/"(L-Af) 



(3.9) 



( Wi\ , f W2\ — 2U"M 

This is an elliptic system for (Wi, W^2)- As mentioned above, if the surface evolves as (13. 6p 
with {Wi, W2) determined by (j3.9p . the coordinate will always be isothermal. 

Remark 3.1 The above system can also he obtained by using 12. 6\) directly. 

Let us conclude this section by deriving the elliptic equations for E and R. Noticing that 
E = Ha ■ Ra = R^ • R^ and Rq, • R/3 = 0, we have 

AE = (Rq • Ra)/3/3 + (R/3 • R/3)Q-a — 2(Rq, • R/3)ci,/3 

= 2(Ra^ • Rq^ — Haa ■ R-/3/3)- (3.10) 

On the other hand, we have 

R„ • AR = 0, n • AR = L + N = 2EH, 

which means that 

AR = 2^ifn. (3.11) 

Remark 3.2 From 13.10\) and the standard elliptic estimate, it is easy to find that E has 
the same regularity as the surface. This fact is noted by S.-S Chern in 0/. Then we can 
gain two more regularities o/R from the regularity of the mean curvature H by using 113. 
Specifically, we will use i3.10\) and \3.11\) to construct our approximate solutions in Section 
5.1. 
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3.2 The velocity equation in the isothermal coordinate 

Assume that u{x, t) is the velocity of the fluid in the isothermal coordinate. Hence, v(u, t) 
u{x{u,t),t), where u is the Lagrangian coordinate. And, we have 

, , ,^ , dR(x(u,t),t) 

u{x{u,t),t) = v{u,t) = - ^ ^ ' ^' ^ 



dt 

_ dn ^ dx^{u,t)d'R ^ dx'^iu,t)d'R ^ 
dt dt dx^ dt dx'^ 

— # 1 — * 2 — ♦ 
= W o X + XjRq, o X + R/3 o X. 

Hence, we have 

x\{u,t) = [-^(u- w) -t^] ox, x1{u,t) = [-^(u- w) -t^] of. 



Consequently, 

dv{u,t) _ du(x(n, 
Ft ~ dt 



du ^ dx (u,t) du ^ dx (u,t) d\i 

— — O X H — — O X H -— r O X 

dt dt dx^ dt dx^ 

^ + 1 [(u - w) . t^] u„ + [(u - w) . t^] U,) O X. 



The above equation can also be derived by Oldroyd's theorem [15j . 
On the other hand, by (j2.2p we have that 



R-a.a = Ln, Rq,,/3 = Mn, R/3,/3 = Nn, 
H = ^{L + N), ArH = ^AH. 

Given that the right-hand side of (j2.15p is coordinate- invariant, the first equation of (|2.15p 
can be reduced to 

^ + ^[(U - W) • + -^[(u - W) • t>/3 

= _2HUn - - + ^ (SuL + 25i2M + S22N) n 

^/E y/E J 

-K^ + 2l^((^-^)' + '"^'))"- ^'-''^ 
Here Sn = Uq, • Rq,, 5i2 = ^(uq, • R^ + • Rq), = • R-/?- So, 

cl _ "^11 q2 _ nl _ '5'l2 c,2 _ 5*22 
1 — -^l — J2 — — , ^2 — 

Qll _ •S'll q12 _ q21 _ 'S'12 ^,22 _ 'S'22 
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3.3 New equivalent system 

Setting u(f , t) = J/^n + Uit^ + U2t^, we infer from that 



/in + 9nti+5i2t2, (3.13) 
- {Up+ — )„ + (t/i^- -_^)t 



=^ /2n + g2it^ + 522t^ (3.14) 
Using (1277]) and ([33]), we find that 







'dt ~ 




ae 




dt 




dn 


1 


'at " 





(^it^ + ^at"). 



Thus, the equation ()3.12p can be rewritten as 



= -2i/nn -^e- + ^(SiiL + 2S^2M + S22N) n 

Ve Ve J 

+ ^^ {Sll,a + 5'i2,/3)t^ + ^^ {Si2,a + 5'22,/3)t^ 

4(^+2l^(^^-^)' + '^''))"- ^'-''^ 

Here we used the fact that E'^{^S^^ + 5^^^) = Su^a + 'S'i2,/3- 

Now we calculate Su^a + 'S'12,^. By (|3.13p and p.l4p . the surface strain rate tensor can 
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be written as 



511 = Uo,-Ra = ^Uia + ^^^- LU"" 

= {Veu,)^ + Mi:^-lu-, (3.16) 

512 = '5*21 = ^ (uq • R/3 + U/3 • Ka) 



= \{^U,), + \{^U,)^-^^^^^^§^-MU-, (3.17) 

= {VeU,), + ^"^^ 7f - NU\ (3.18) 
2/E 

Since the incompressible condition Vr • u = can also be written as + ^22 = 0, we get 
by (|3lB and (|3T8D that 

{VeUi)c, + {VeU2)p - = 0. (3.19) 
We get by ([331) and + 522 = that 

5ll,a + 5i2,/3 = Siia — 2r\iSii + Si2(S —^22^li ~ ^12^12 
= Siia + 5i2/3. 

Similarly, 

5l2,a + 822,(3 = 5i2o + 522/3- 

Using (f3J6]l - (t3J9]l . we find that 

5n« + 5,2, = ((^/^[/l). + ^^%^|^-L[/" 
^ 2v ii/ 



+ 



l(^/lf/,), + i(v/lf/.)„ - - MU" 



^(:^+lf(V£f/,)„ + (^/l(/,), 



2VE 'o. V 2\/^ 

2 V 

^£;<,c/2 + ^/3C/i 



2/E 



+ MU' 
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Thus, we obtain from (|3.15|) the evolution equation for Ui: 



L[7" + EHW' 



dt E^/E y/E 

2eo / Efi{J2 — Eq[J\ 



E^\ 2VE 



2v^V /a Ey/E\ y/E 

^ 'U2Wi^ + iUi-Wi)Uia + iU2-W2 



'E 

+i(2C/2C/"M + UiU^L) - — L=((C/i - Wi)U2E^ - U^E^ . (3.20) 

A similar evolution equation for U2 can also be obtained, although we omit the details here. 
The evolution equation for the normal velocity U"' is 



^AH -2HU-^ ((2C/i - Wi)U^ + {2U2 - W2)U^ 



at/" 

"aT ^ 2^"" s/e 

E\ E 

+-^{LEpU2 - ME^U2 - MEpUi + NE^Ui) 
E'^y E 

2enC/" ,^9 . „, ,2 



(L^ + 2M^ + Af' 



^ ^(L- A)2 + 4M2) - l(LC/f + 2MC/1C/2 + AC/|). (3.21) 



4:E'^V ' I E 

Due to (j2.9p . the evolution equation for the mean curvature H is 

1 , VFi , , C/" ..2 , o,.r2 , ..2 



We denote 



EpU2-E^Ui U^.^^ ^. ^ _ 1,^^„,2 



F, = ^^^^ + -(A-L), F2 = -(^") 



F4 = -^U2Wi^ + i(2M[/2[/" + Lt/if/") - — ^ (([/i - VFi)C/2B/3 - U^E^ , 
V -C/ -C/ 2E\J E ^ ' 



F5 = -^UiW2a + ^{2MUiU^ + AC/2C/") - — ^ ((C/2 - ^2)f/lSa " U^Ep , , 

V-B E 2E\/E V 



° (L£;;3C/2 - MBc.C/'2 - MEpUi + NEaUi) ^{L^ + 2M^ + N^) 



B2/E B^/^E 

/ \ 1 

(L - A)2 + ) - -(LiJi^ + 2MU1U2 + AC/|), 



4B2V' ' J E 

Tjn 

Fr = ^{L^ + 2M^ + N^), 
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and Ui = Wi = ^ for i = 1, 2. 

Prom (j3.19p - (|3.22p . we obtain the following equivalent full system: 



C/"n + VFit^ + Wat^, (3.23) 



dR 

'dt 

-{Ui - Wi)Uia - {U2 - W2)Uip + F4, (3.24) 
-(til - Wi)U2a - {U2 - W2)U2i3 + -^5, (3.25) 

^ = -^^^ - - ((2ni - wi)U: + {2u2 - W2)U^[ 

+ -^{LUia + MUip + MU2a + NU2f3) + Fe, (3.26) 

f)ZJ 1 

— = —AU'' + wiHa + W2H^ + F7, (3.27) 
(\/;BC/i)« + {VeU2)p - 2EHU'' = 0, (3.28) 
where iWi^ W2) is determined by the elliptic system p.9p . 

Remark 3.3 Actually, I3.21\ j is induced by i3.23\) . However, in order to perform a suitable 
energy estimate, we add it to the full system. 

3.4 The equation of the pressure 

Using the incompressible condition {\fEUi)a + {VEU2)f3 = 2EHU"^, we find that 

{VEUu)a + {VEU2t)f3 - 2EHUr = U''{2EH)t - (:^U,) ^ - (:^U, 
We denote the left-hand side of (j3.15p by 



2/E ^« 



Noting that Et = 2\/!EAoi, then we have by (f3TT5]) that 



E 

-Kio, - + 2^HK'' - {SiiL + 2S12M + 522 A^) 

+H (^AH + ^{{L-Nf + 4M2 
= U''{2EH)t - (AoiC/i), - {AMp- 



14 



By a direct computation, we obtain 

U^{2EH)t + (Ki - AoiC/i)a + {K2 - AoiU2)p - 2VeHK^ 

= (U,U,^ + U,U^, - ^{2MU2 - (L - N)U,) + ^i^^^^J^I^^ 



' ' ' 2E Jp 

-4HVe{UiU'^ + f/sC/^) - 2H{LUf + 2Mf/2C/2 + NUl) - {U^f - {UJ^f 

+ (E!)!(l2 + + - 4i?2^2) def 



And using the incompressible condition again, we get 



5*110 + 'S'21/3 \ / 'S'i2a + S22j3 \ 

E )a \ E ) H 

-^(5n, + S21P) - ^iSi2a + S22P) + ^Ai2EHUn 
+4 ({^)fsU2 - {^)JJi - LUA - I ({VEUU2 + {VE)fsUi + MU' . 



Consequently, we obtain 



-An + AEH^U = h(aH + ^({L- Nf + 4Af 
V 2E 



+Gi - 2eoG2 - {SuL + 25i2M + S22N) =^ g. (3.29) 
E 

Remark 3.4 In (MW> G is a polynomial function of{d^E, d^U'\ d^Ui,d^U2, d^L, d^M, d^N, d^H), 
where < < 3 anc/ < |/| < 2. 

Remark 3.5 It is reasonable that there is no term involving w in Q. Actually, by differen- 
tiating the equation a° • v^^ = 0, and reformulating the resulting equation in the isothermal 
coordinate, we can also derive the equation of the pressure. 

4 The linearized system 

In this section, we study the well-posedness of the linearized system of (|3.23p - p.28p . More 
precisely, we will consider the linear system 

f ^ = -^AiVEUi) + G„ i = l,2, 

9U^ = -^AH + B^-VU'' + G3, (4-1) 

together with the initial condition 

(f/i, U2, f/", H)\t=o = (f/i°, Ul U^, i^O). (4.2) 
Throughout this paper, we assume that x = {xi,X2) G T^. 
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Theorem 4.1 Let s = 2(k + 1) for some integer k > 2, and let T > 0. Suppose that E G 
C([0,T];F'*+i(T2))nCH[0,T];i7^-i(T2)) and E > cq for some cq > 0. We also assume that 
{UlUlU^,H^) G i7-i(T2), (Gi,G2) G L2(0,r;//-2(T2)),(G3,G4) G L2(0, T; F-^T^)), 
and iB^,B'^) G L2(0, T; i?''^^(T2)). T/ien i/iere existe a unique solution {Ui,U2,U'^, H) on 
[0,T] to the linear system OTTP-OT^ such that 



{Ui,U2) G C{[0,T]-H'-\T^))nL\0,T-H'{T^)), 
G C{[0,T];H'~\T'^)). 

Moreover, for any given e > 0, it holds that 

Es{t) < C{\\E\\loohs-i)[Es{0) + ^ J^e{\\E\\Hs+^,\\Et\\Hs-^){l + \\B\\Hs)Es{T)dT 

+ f \\{G^,G2){T)fH.-.dT + e f\\{G3,G,){T)\\l.^,dT], (4.3) 
Jo Jo 

where B = {B^,B'^), is an increasing function, and Eg{t) is defined by 

E,{t)tf\\{U^,U2m\\j,..,+ [' \m,U2){r)\\hdr + \\{U\Hm\\h-^- 



Proof. The existence of (Ui, U2) is ensured by the classical parabolic theory, whereas (C/", H) 
can be obtained by the duality method, see [l] for example. Here we only present the proof 
of the energy estimate. For this purpose, let us introduce the energy functional £ defined by 

with £^ and £'^ given, respectively, by 

£'^M\\A^-^U^\\l, + \\A^-'U2\\h, 

£' ||a(1a)'=C/«||2, + \\A{^A)''H\\l,, A = (-A)i 
Step 1. Estimate of £^ 

Taking the derivative to £^ with respect to t, we obtain 

—£' = {A'-'Ui,A'~'dtUi) + {A'-'U2,A'-'dtU2). 
By using the first equation of (j4.ip . we get that 

{A'~^Ui,A'-^dtUi) = {A'-'^Ui,A'-'^^^A{^Ui)) + {A'-'^UuA'-^Gi). 

Ey E 

By the Cauchy-Schwartz inequality, we have 

{A'-'Ui,A'-'Gi) < ||C/i||j^.||Gi||^,.-2, 

and we write 

(A-i[/i,A-i-^A(^/^[/i)) 

= (A-if/i, [A-\ ^^^] Ui) + (A-i[/i, -^A(V^A^-iC/i)). 
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By integration by parts and based on Lemmas 16.11 and 16.21 the second term of the right-hand 
side is bounded by 



-c||A^C/i||i.+C7||Vi?||i.o||A^-iC/i||i.+C||Vi?||Loo||A-i[/i||i.||A^C/i|U2 

<-l\mh+mE\\H^^^)\\ui\\is-.. 



Here the constant c > depends only on cq. And when Lemma 16.31 and Lemma 16.21 are both 
used, the first term is bounded by 



since we can write 



Summing up the above estimates yields that 

{A^-'Ui,A'-'dtUi) < -c\\UifHs+T{\\E\\Hs+i)\\Ui\\j,.., + \\G,\\l..„ 
for some c > 0. Similarly, we have 

{A'-'U2,A'-%U2) < -Cimh + H\\E\\Hs+^)\\U2fHs-^ + 

Hence, we obtain 

p' + c\\{UuU2)\\jjs<Tm\Hs.^mUi,U2)\\l^^^ (4.4) 
Step 2. Estimate of S"^ 

Take the derivative to with respect to t to obtain 

= {A{^A)%U\A{^A)'^U-) + {A{^A)%H,A{^A)'H) 
+ {A[du{^A)'']U-,A{^A)'U^) + {A[dt,{^^)']H,A{^A)'H) 
# / + // + /// + IV. 
Based on the last two equations of (|4.1|) . we get 
/ + // 

= (A(l A)^(i?i • Vf/"), A(i A)'=C/") + (A(l A)^G3, A(i A)*^i7") 
+{A{^A)'{B^ . Vi?), A(lA)^if) + (A(-iA)^G4, A(lA)^if) 

=^/l +/2 + //1 +//2. 

Here we use the following fact: 

{A{^A)H-^A)H,A{^A)'U-) + (A(1a)^-(^A)C/", A(lA)^if) = 0. 
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Using Lemma 16.61 and Lemma 16.41 we get 

\iii\ + \iv\ < H\\{E.Et)\\Hs-.){w\\is-. + ml..;), 

and by Lemma 16.41 

1/2I + 1//2I < mEWus^^Mu'^wis-. + +^11(^3,^4)11^.-1. 

To estimate /i, we write 

E EE 
+{B' • VA(^A)'=C/", A(lA)'=f/") = 111 + I12. 

We have by Lemma 16.41 that 

\Ii2\<H\\E\\Hs-^)\\VB^\\L^\\Ums-^. 

We further write 

111 = (A[(-A)^Sl] • VC/"+ [A,^!] • V(-A)'=C/",A(-A)^[/'^) 
E E E 

+(^1 • A [(i A)^ V] A(iA)'^t/"), 
which along with Lemma 16.41 and Lemma 16. 6116. 71 imphes that 

\hl\<F{\\E\\Hs)\\B^\\Hs-.WfH..,. 

On the basis of the above estimates, we obtain 

^ d „2 / Tf\\Tp\\__ 1177 11 _\(^ , IIDII .\{\\TTn\\2 , n ,^112 



£^ < J-.(||i?k.,||ii;tib.-i)(i + ||i?||H-0(r"llH-i + ll^ll?/-i) 

+e||(G3,G4)|||. 



Js-l- 



Step 3. LP' estimate 

Taking the energy estimate for Ui{i = 1,2), we obtain 



< C(||i?||i^)||-^Vi?||L^||V(^/:Et/0llLH|t/.||L2 + \\G,\\lA\U^\\l^ 

< H\\EM\\U^\\h + WGiWh + ll|v(^/^C/.)||i.. 
Taking the energy estimate for {U^,H), we get 

{EdtU'',U'') + {EdtH,H) 

= {EB^ ■ VC/", U") + {EG3, U") + {EB^ ■ VH, H) + i^), 
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from which, we infer that 
1 d 



< J-,(||(i?,i?0lk3)(l + ||B||^3)(||C/"||i, + ||/?||i.)+^ll(G3,G4)||i.. 
Thus, we obtain 

po < n\\mH-^mui,u2)\\i.+M\mEt)\\Hs){i + \\B\\^^^^^^^ 

+ \\iGi,G2)\\l.+e\\{G3,G,)\\l,. (4.6) 

Here So =^ ?72)||i. + WVEU^I, + \\VEH\\1,. 

Now we are in position to complete the proof. Taken together (j4.4p - (j4.6p yields that 



di 



{8 + r]£o) + c\\iUuU2WHs < \\Et\\Hs-i){l + \\B\\hs){£ + ^£o) 

+77(||(Gi,G2)|||.-2+£||(G3,G4)|||.-i), 

which implies (I4.3P by taking rj to be bigger than C{\\E\\i^aaf^s-i), since we have by Lemma 
6.51 and an interpolation argument that 



>c\\iU",H)\\l.^, -Cm\Hs-^)m^H)\\l,. 
This completes the proof of Theorem 14.11 □ 

5 Nonlinear system 

This section is devoted to solving the nonlinear system (I3.23p - (l3.28p . 
5.1 Iteration scheme 

We will construct the solution (Il,Ui,U2,U"' , H) by the iteration method. First of all, we 
take 

And, {Wi^\W2^^) are determined by solving the following elliptic system: 

^{0) \ / „^(0) X 2!7"'(0)M(0) 



E(0) 



( wr \ ( wr \ 

The pressure is given by 

-AnW +4£;W(F(o))2n(o) = 

with determined by (R(°), C/f \ C/2^°\ f/"'(°)) (see I^M))- 
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Assume that (C/f \ ;7^^\ P^f \ \ RW) has been constructed. We de- 



note 



LW=RW-nW, MW=R(f).nW, iVW=Rg).nW, 



Then we construct {u[^'^^\U2~^^\u^''^^~^^\ H^^~^^^) by solving the foUowing hnear system: 



■/3 I 



at 



eo 
eo 



1 A^./F??^rr(^+l)^ I x-W 



Qljn,(i+1) 



w 



(5.1) 



where u 



c/f) (^) 

J 7/;: ' 



for i = 1,2 and 



2eo 



pie) , 



3/3 



2eo 



eo 



1 Jiffl 

2\/:ew 



3/3 



+ 



2eo 



6 ' 



with = I,-- - ,7) are given in Section 3.3 where {Ui,U2,U'' ,Wi,W2, E, L, M, N) are 

replaced by (C/f ^ , C/^^^ , C/'^'^^) , VF^ ^ , , E^^) , L(^) , M^^) , 7V(^) ) . Let R'^+^ be given by 



R 



Rl 



t=o 



Rn. 



And, R^^"^"*^) is determined by solving 



AR(^+i) - r(^+i) = 2i?WR(f+i) X r('+^) - r(^+i). 
Then we construct the surface R^^"^^) by solving the following elliptic equation: 
AR(^+i) - R(^+i) = 2i?(^)R(f+i) X R(f+^) - r(^+i). 
Next we define by solving 

A^(^H-i) - 2ii;(^+i) = 2(Rl^;i).R(7^)-R(f+i).R(7^)) 



(5.2) 



(5.3) 



(5.4) 



fR(f+i) • R(f+i) + R 



■/3 



'/3 



(5.5) 
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And, {Wi^^^\W2^^^^) is determined by solving 



(5.6) 



Finally, we define the pressure n^^"*"^) by solving 

- An(^+i) +4i?W(/?W)2n(^+i) = gw, (5.7) 

with determined by ([/f \ C/f \ C/'^'W, ^W, L^, M^, iV^), see (fCT]) . 

Remark 5.1 //R'-'^"'^^ is directly defined by i5.2\) . then we can only obtain the H^~^ regular- 
ity o/R'^'"'"^). However, we need the H^^^ regularity ofR^^^^^ to close the energy estimates. 
Motivated by 13.11\) . we determine R^'"*"^^ by using / f5.^) -( f5^ so that the H^'^^ regularity of 
can be obtained by the elliptic estimates. 

5.2 Nonlinear estimates 

Before presenting the estimates, let us make the following assumptions on the step-^*'' ap- 
proximate solutions (R(^) , t/f ^ , /7f ^ , /7'^'(^) , E^^^ , H^^^ ) : 

sW(x,t) > Co > 0, for any (t,x) G [0,r] X T^, (5.8) 



(5.9) 



/ (i?(^)(x,t))2(fx > ci > 0, for any t G [0,r], 
|R(f)(x,t) X R^^^(x,t)| > Co, for any (t, x) G [0, T] X (5.10) 

||(^f\C/f ||ioo(o,T;H-i) + ll(f/f\C^f )llL2(0,r;H^) < Ci, (5.11) 
||(^"'W,FW)||io.(o,T;H^-l)<C2, (5.12) 

l|R-^^^llc*([o,T];//''+i-2«) + ll-E'^^^llc"{[o,T];H''+i-2') < C's, i = 0, 1, (5.13) 

\\^^^^\\l°°{Q,T;H''-^) + \\E^^^\\l^(0,T;H=-^) < C^. (5-14) 

Here T > 0, s = 2{k + l),k > 2, and Ci,C2,C3, and C4 are some fixed constants to be 
determined in Section 5.3. Note that the assumptions (|5.8p - (j5.10p and (|5.13p - (j5.14p are made 
so that we can use Theorem 14.11 at each step of the iterations, and the assumptions (I5.1ip - 
(j5.12p are determined by the energy estimates for the linearized system. 

In what follows, we denote C by an increasing function, which may be different from line 
to line. From the definition, it is easy to see that 

||(LW,MW,ivW)||i.o(o,T;H-i) <C(C3). (5.15) 

Using Lemma 16.81 and Lemma 16. H we find that 

ll^^f^'^llL-(0,T;H^) + l|Wi'^'^llL-(0,T;H=) < ^(^2,^3). (5.16) 

From ()5.15p . ()5.16p . and Lemmas I6.m6.21 we infer that for z = 1, • • • ,7, 

\\Fi''^\\L^{0,T;H'-^) < C(Ci,C2,C3). (5.17) 
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Thanks to Remark 13.41 we get by using Lemma |6. m6.2l that 

ll^^^^llL°°{0,r;_f/''-3) < C{Ci,C2,C3). 

Thus, we infer from Lemma 16.91 that 

l|n('"'^)||L^(0,T;H-l) <C(Ci,C2,C73). (5.18) 

Using (IST^ - (I^3HD . we obtain 
Proposition 5.2 The nonlinear terms F^\F2 \ and f'^'^ satisfy 

\\^f\\L°°{0,T;H'-'^) < C(Ci, C2, C3), i= 1,2, 
II-^P||l2(o,T;H''-1) < C(Ci,C2,C3). 

In order to prove the convergence of the iteration scheme, we need to establish some 
difference estimates in the lower-order Sobolev spaces. For this, we set 

4 = uf^ - Ut'\i = 1,2), 6fj. = C/-'W - 4 = - F(^-i), 

4 = rW - r(^-i) , 4 = - i^(^-i) . 

First of all, we have 

||(LW,MW,ivW) - (L(^-^),M(^"^),Ar(^-^))||^.-3 < C(C3)||4k=-i, 

livt^(m) _ < c(C3)(||4k»- + + II4"IIh=-3), 

which imply that for i = 1, • • • ,7, 

lli^f -i^^'"'^llH»-3 < C(Ci,C2,C3)(||4k^-i + ll4llH»-i 

Similarly, we can obtain 

||gW_g(^-i)||^,„3 < C(Ci,C2,C73)(||4lb^- + ll4llH»-i 

Hence, we infer from Lemma 16.91 that 

lin(^+i)_nW|l^^_3 < ciCi,C2,Cs){\\si\\Hs-i + \\6U\Hs-i 

+ ¥ui\\h=--^ + ¥u2\\h=--^ + ¥u"\\h=-'-' + ¥H\\H=~'i)- 

From the above estimates, we can deduce 
Proposition 5.3 For i = 1,2, it holds that 

< c(Ci,C2,C3)(||4k-i + ll4llH-i 

-J't'^WHs-s < C(Ci,C2,C3)(||4llH-i + ll4lb-^ 
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5.3 Proof of the main result 

To simplify the analysis, we will first prove the well-posedness of the system by assuming 
that the surface can be globally parameterized by the isothermal coordinates. In Section 5.5, 
we will indicate how to extend this result to a general closed surface, and thereby conclude 
the proof of Theorem II. 1[ 

Theorem 5.4 Let s = 2{k + 1) for some integer k > 2. Assume that U2, Uq) G 
//*~^(T^), and the initial surface Rq G H^'^^ . Moreover, the coefficient of the first fun- 
damental form Eq and the mean curvature Hq satisfy 

{y%Uf)a + {y%U^)p - 2EoHoU^ = 0, 
Eoix) > 2co, / Hl{x)dx > 2ci, 

for some cq > 0, ci > 0. Then there exists T > such that the nonlinear system i3.23\ }- [3.28\} 
has a unique solution (R,Ui,U2,U"' , H) on [0,T] satisfying 

(C/i, U2) e C{[o, T]-w-^) n L2(0, T; H'), 

RG C([0,r];i7^+i), G C{[{),T\-H'-^). 

Remark 5.5 We have chosen the isothermal coordinate for the initial surface. Hence, the 
conditions 

Eoix) > 2co, Jrj.2 Hl{x)dx > 2ci 

are naturally satisfied for any smooth closed surface. 

Proof. We split the proof into two steps. 
Step 1. Uniform estimates 

Let us assume that (RW, C/f \ C/2^^\ F^) satisfies ^^-<^J^. We will show 

that (R(^+i),C/f+^\c/2^^+^\c/"'(^+i),£;(^+i),i7('^+i)) also satisfies the same estimates. 
We denote 

4^)(t)'= + /^l(f/!^f/f 

JO 

Then we infer from Theorem 14. II and Proposition 15.21 that 

E'i+'-Xt) < C{C^)[Eo+C,{Ci,C2,Cs)l\l + \\{ui'\u^)\\Hs)Ei'+'Hr)dT 
+C{Ci,C2,C3){t + e)y 

Here Eq =^ \\{Uf, U^, U^f , Then we get by Gronwah's inequality that 

E^^-^^\t) < (C(C4)i^o + C(Ci,C2,C3,C74)(t + e))exp(C(C7i,C2,C3,C4)t). 
Taking T and e small enough yields that 

4^+1) (t) < 2C{Ci)Eo for t G [0,r]. 
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This means that if we take Ci = C2 = 2C{C4)Eo, (C/f C/^ ?7"'(^+i), satisfies 

dsn])- dsn]). 

Due to (|5.2|) . we find that 

Hence, by taking T to be small enough if necessary, we get 

\\^^^~^^^\l°°{o,T;H'>-^) < 2||Ro||/i-s-i. (5.19) 

We also have by (|5.2p that 

||atR(^+i)||ioo(o,T;/f-2) < C{C2,Ci). (5.20) 
We get by the elliptic estimate that 

\\^^^~^^'^\\l°°{0,T;H'>) < C(||Ro||//s-i, C2), 

which along with (j5.19p and (15. 4p implies that 

\\^^^~^^'^ \\l°°{0,T;H''+'>-) < C(||Ro||_H-s, C2). 

Hence, by (j5.5p and the elliptic estimate, 

l|-^^^'''^^llL°°(0,T;_ff=+i) < C(||Ro||//», (^2)- 

Taking the derivative to (j5.4p and (|5.5p with respect to time, we get by (j5.20p that 

||9tR(^+l)||i.o(o,T;H-i) + ||5t^('+^)||ico(o,T;i^-i) < C(||Rolk»,C2,C4). (5.21) 

Hence, taking C3 = C(||Ro||h-, C2, C4), we see that (R'-^+^\ E^^+'^^) satisfies (l5l^ . 
Now taking C4 = 2(||Ro + ||^o||h=-i) < C\\Ko\\h->, it follows from (fOT]) that 

||(R(^+l),i?(^+l))||i^(0,T;H-i) < ||Ro||h-i + ||i?0||//-i +C(||Ro||//.,C2,C4)r, 

which implies that (R.^^^^\ E^^~^^^) satisfies (j5.14p when T is taken to be small enough. 
Similarly, we can show that (R(^+i), ^(^+1), if(^+i)) also satisfies (ISlSjl - dSTOD . 

In conclusion, we prove that there exists a T > depending only on \\{Ui, U2, Uq)\\hs-i 
and ||Ro||h»+i such that ^-^M hold for iR(^+^\ui^-^^\u^^^^\u'<^+^\ E^^+^\ H(^+^')). 

Step 2. Existence and uniqueness 

It suffices to show that the approximate solution sequence is a Cauchy sequence. For this 
purpose, we set 

5^+1 = [7.(^+1) -[/W(i = 1,2), 5^^+i = [/"'(^+i)-[/-'W, 4+i=F(^+i)-/fW, 
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Then {6^^^ , 6^^^^ , dp'J' , 6^j^^) satisfies the following system: 



eo 



in 



( 95^ 

at 

85'+^ 
dt 

dt — 2EW 

dt ~ 2E^i) '^Ojjn -f- W- OiO^ -t- , 

L {5'u';\sX,sp^',s'^%=o = (0,0,0,0), 



(5.22) 



where 



6t: 



(i-i) 



1 



1 



1 



1 



A(/^(^C/f)), 



5f: 



F, 



+ eo ==A(\/^C/f ) - eo ; 

_ ^Ai^W + ^-A//W - {25l - <5l.)a,C/"'W, 



A(Vi^(^-i)C/f ), 



7 "7 ' 
From Proposition 15.31 it is easy to see that 

¥Ff\\Hs-. < C{\\6i\\H^-. + \\6%\\hs-i + \\{6fj,,Sfj^,Sf;^,Sjj)\\Hs-s), i = 1,2, 
||5J;('^||h=-3 < C(||4||^.-. + 11411^.-1 + ||(5^„4)||h^-. + \\{Sfj.,6U\H^-.), 
Revisiting the proof of Theorem 14. H we can obtain 

D{+\t) < Ce f D[^\T)dT + cj2 f ¥Ff{r)\\l.^,dr 
Jo Jo 

+Cej2 f ¥F\\r)\\l..,dT 
i=3 -^0 

< Ce [ D{+\T)dT + C{t + e) sup D^{t), 

Jo re[0,t] 

where D^{t) is defined by 

D'{t) = Di{t) + \\5imH^-^ + \\6'E\\H'-^, 

with Diit) = ||(<5f,,,5^,,5f;.,5f,)(t)||H^-3 +/o*||(5^^,<5f,J(7 



3,4. 



(5.23) 



ydr. 



On the other hand, we revisit the proof of Step 1 to find that 

\\&\hs~^ + K^'Whs-^ < C,D[{t) + Ct{\\5i\\ns-. + \\5'e\\hs 



(5.24) 



For some small (5 > depending only on Co, with e and T taken to be small enough, it follows 
from (fOSD and dOijl that 

sup {D[+\t) + 5Di+\t)) < I sup iDiit) + 6Diit)), 



te[o,T] 



te[o,T] 
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with = ||<5^^||h=-i + \\S^e^\\hs-i. This imphes that 

is a Cauchy sequence. More precisely, there exists the hmit {'R,Ui,U2,U'^ , E, H, 11,11,11) 
such that 

rW^R, E^^^ ^ E in L°°{0,T;H'-^); 

Uf^ ^Ui, u!i^ ^U2 in L°°{0,T;H'~^)nL'^{0,T;H'-^y, 

rW^R in L°°{0,T;H'-^), R^ ^ R in {0,T; H'~^). 

With the above information, it is easy to prove that (R.,Ui,U2,U"' , E, H, 11,11,11) satisfies 
the system (j5.ip - (|5.7p without the index i. In particular, we have 



^ = ^R„xR, + ^R„ + ^R„ (5.25) 



AR - R = IHKa X R/3 - R, (5.26) 

AR - R = 2HB,a X R^ - R, (5.27) 

AE-2E = 2{-Ra^ ■ R„^ - R«^ • R^^) - (R^ • R« + R^ • K/s), (5.28) 

f - - - 7i"^ - w^'' - - <^-^»> 

L = Raa-n, M = R^^-n, iV = R^^ • n, n = (5.30) 

\ti.a X il^l 



And, iWi,W2) satisfies 



( f mS\ - f m.\ - u"iL-N) 



(5.31) 



It remains to show that the solution of the limit system is a solution of the original system. 
For this purpose, it suffices to prove the following relations: 

R = R = R, E = Rq, • Rq, = R^ • R^, R^ • R^ = 0, H = — . (5.32) 

And the incompressible condition follows easily from (|5.32|) . As the proof is very complicated, 
it will be given in the following subsection. 

5.4 Consistency with the original system 

This subsection is devoted to proving (I5.32p . Let us introduce some notations: 

ail = Ra • Raj ai2 = Ro • R/3, 122 = R/3 • R/3, 
Oil = Rq • Ra, ai2 = Rq • R/3, a22 = R/3 • R/3, 

w — — Rq X R^ H — ^Ra H — i="-fSi 
E \/E VE 



\^ = {wi,W2), Wi = ^^{i = 1,2). 

V E 
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We set 

(5_R = R-R, (5^ = R-R, 5E = E-an, 6^ = E - a22, = an -0.22- 

In what follows, we denote by 5 some operator bounded in iJ*''(T^)(0 < A; < 1), which may 
be different from line to line. For example, 

\\'S{u,v)\\Hk < C{\\u\\Hk + \\v\\Hk). 

We get by using (fOSD and (fOTT) that 

(an - 022)4 = w„ • Ra - W/3 • R/3 

= w„ • R„ - w;3 • R;3 + 5(V(R - R)) 

C/"(L-iV), „ , ,1, 

= (2|Ra X R^l - On - 022) + -("^1" + W2i3)[aii - 022) 

+wi{aii - 022)0 + W2{aii - 022)/? + [wia - Wii3)ai2 + 5(V(5h,, VJ^). 

Similarly, we have 

(012)4 = • R/3 + • R„ 

= w„-R^ + W;3-R« + 5'(V(R-R)) 

MV^ 1 
= ^ (qii + 022 - 2|Rq X R^jI) + -{wip - W2a){aii - 022) 

+Wiai2a +W2ai2l3 + ("i^la + W^2/3)ai2 +^{^5r,V5~). 



Noting 



(an + 022)^ - 4(ana22 - af 



flu + 022 — 2|Rq X R^l = ^ ^ ""^ ^'^ 



On + 022 + 2|Ra X R_a| 
2 

12 



(an - 022)^ - 4a2 



Oil + 022 + 2|Ra X R^ 



we find that 

dtda + ^-V5a = d{ySR,V6j^,V^6R,V^6j^,6a,ai2), (5.33) 
dtai2 + ^ ■'^ai2 = d{ySR,V6j^yHR,VH^,6a,ai2)- (5.34) 

By ([Oe]) and (fOTl) . we have 

(A - 1)(R - R) = 2i?(R„ X % - Rc, X R^) = 5'(V<5^). (5.35) 

And by ([CT]) and dS^H]), 

(A -1){E- an) = -2(Ra • {2HRa x R^))^ + ^ _ 022 - 2(R„ • (R - R))«, 

which implies that 

{A-2)6l, = d{V6R,V6^,6a}. (5.36) 
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Similarly, we have 

(A-2)(E-a22) = -2(R;3-(2i?R„ xR;3))^ + 022 -aii-2(R^-(R-R))^ 

= d{ySR,V6^,6a), (5.37) 

(A - l)ai2 = 2H{{Ka X R„^) • R^ - (R« x K^p) ■ R^) 

+ 2H{{KalS X R;3) • Ra " (Ra/3 X Kfs) ■ R^) 
2 

= 5^5(V%,V'=(^^,Si2). (5.38) 

k=0 

The following facts will be used frequently: 



|R« xR^I - ^ = y^det(a,^)-S = i?(5^,j|,ai2), (5.39) 

r}i, r?2, ri,, -r^^ = ^ + ;?(5i 4, vJi vj|, ais) , (5.4o) 

r}2, r^i, ri2, -r?i = ^ + 5(<5i, <5|, v^i v4, 012) . (5.4i) 

{AR)^ = i2HRaXKp)^ + d{ySRy6j^,V^6R,V^6j^), j = a,/3. (5.42) 



Indeed, we have 

■n2 1 / ai2 o , 0.11 r. \ 

= |R.xR,| (-^^^"^^ + ^^""^^) 

= ^+5(<5i<5iV<5|,«i2), 
and the others can be deduced similarly. For the last fact, we have by ()5.27p that 

AR = 2HRa X R;3 - R + R 

= 2HKa xKi3 + 2H(R.a X R^ - Rc, X R^) - R + R, (5.43) 

thus (|5.42p follows easily. 

To proceed, we also need the following lemma. 

Lemma 5.6 For ^ = a, /3, it holds that 

2 1 

{2EH -L-N)^ = J2Y1 ^l^"'*^^' ^'■'^5' "^'^E, V^ai2) , (5.44) 

fc=0 j=0 

1 /?2 I p2 2 

LN -M' = ^[-AE + + E H'^'SIVHI, V'^au) , (5.45) 

k=0 

N^ = Mp + HEa + 6l,V6l,V5l, au) , (5.46) 

L^ = M^ + HEp + 5(<5i 5l,V5l, V^, 012) . (5.47) 
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Proof. First of all, a direct calculation gives 

E 



2EH -{L + N) 



AR • (R„ X R^ 



1 



Rry X Rf 



r AR • (R„ X R^) 



AR 



iJ/Ro, X R^ Rq, X R^ 



|R«X%|2 |RaXR^|/' 

which implies ()5.44p . From the Gauss equation, we infer that 



LN -M^ 





"11/3 
2 

0220 
2 





«12a 


12 y 




«11;3 
2 


2 


an 


ai2 


Ol2 


122 



1 . 



1 . 



1 



1 



2^22/3 



On 

ai2 



0.22 



which implies (I5.45p . The Codazzi equation 6^^ 



7 — "a7„ 



implies that 



dh 



11 



5j;2 



^12^11 — r?26i2 — r|2feii — r?2^2i 



- Fai^ll - r2i6i2 - ^\lhl2 - rfi622, 

rj2^2i — ri2^22 — r22^ii — r22^21 



~ ^21^21 — r2ife22 — r21^12 — J- 21022, 



r^i62 



where bu = L, bu = 621 = M, and 622 = N. Then ([536])- ([5311) follow easily from I^5l0\i 
and (I5.4ip . The proof is completed. □ 
In the following, we calculate R — R. By (j5.25p and (j5.26p . we have 



(A - 1)(R - K)t = Aw - {2HTia x R^)*. 
Direct calculations yield that 



(5.48) 



jjn 

Aw = A( — Ra X R/3 + WiKa + lf2R/5 

E 



jjn 



jjn 



Jjn 



A( — )R« X R^ + -_(AR)« xR^ + — R^ x (AR)^ 



E 



E 



+Au;iR„ + Au;2R/3 + ?^i(AR)„ + w;2(AR)^ + 2— (R„^ - R^^) x R„^ 



+2 



+2 



Jjn 

~e' 

Jjn 



Raa X R^ + 2 I — 

a \ h/ 

R„fl X R^ + 2 



Rq X Rc/3 

Rq, X R^^ 



E )fi "'^ ^ \E 

+2'U;iaRaa + 2'U;2/3R/3/3 + 2{wx^ + W2a)R«, 



(5.49) 
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and by ([225]), 



2i7(R„ X R^) 

/ A/7" [/" \ 

-- (-^ + 2wiH^ + 2w2H,3 + + + iV2)jR„ x R^ 

+2Hi^—K^ X R^ + wiR^ + W2R/3j^ X R/3 

+2i7Ra X (^-^R-Q X R^ + WiKa + 1^2^,13 

/ A/7" /7" \ 

: (^-^ + 2u;i//« + 2w2H^ + -^(L2 + 2Af2 + Ar2)jR^ x R^ 

+2H [ —Ka X K/S + WiKa + W2'Ri3 ] X R^ 

/[/" \ 
+2/fR„ X — R^ X R/3 + t/;iR„ + t/^aR/? + 5(V5k, V(5j|) 



/ A/7" f/" \ 
= (-^ + -^{L^ + 2M2 + iV2)jR^ X R^ + w;i(2ifR« x R^)„ 

+Zi;2(2i^RQ X R^)_g + 2H{wia + t(^2/3)Ra X R^ 

/ /7" \ /7" 
+2i/(^— j^(R„ X R/3) X R^ + 2i/— (R„ x R^)^ x R^ 

/ /7" \ /7" 
+2i7(^— j^R„ X (Ra X R^) + 2i/— R« x (R« x R^)^ + 5(V(5ij, V(5jj) .(5.50) 

From (iO0]) - (fOT]) . it follows that 

Raa = TJ^^Rq, + T\-J\.ii + Ln 

= ^Ra-^R/3 + Ln + 5(5^,4,V5iv4,ai2), (5.51) 
R„^ = rlsRa + r?2R/3 + Mn 

= ^Ra + ^R/3 + Mn + ^(5l„5|,V4,V5|,ai2), (5.52) 

R/3/3 = r22Ra + ^22^13 + 

= -^Ra + ^R/3 + iVn + ^(5i4,V5iv4,ai2). (5.53) 

And by ([5^ . it follows that 
(AR)„ X R^ 

= 2i/„(R„ X R^) X R;3 + 2F(R„ x R^)„ x R^ + ^{V5r, VJ^j, V^^^,, V^^^) 

= -2EH^Yi^ + 2H{Yi^xYip)^xB.p + :S{y5R,V5^,VHR,VH^,6E,ai2), (5.54) 
R„ X (AR);3 

= 2Hpna, X (R„ X R^) + 2//R„ x (R„ x R^)^ + ^{V5r, V6j^, V^6r, V^5j^) 

= -2EHp^p + 2HR^ X (R„ x lip)^ + ^(VJ^j, V5^, VHr, VH^, 5e, a^) ■ (5.55) 
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Summing up (j5.48p - (j5.55|) . we obtain 

(A - 1)(R - R)t = Aw - {2HR.a x R^)t 
= Al^—jEn - 2C/"(i7„R^ + HpTi^) + AwiK^ + Aw2'Rp 

+2(— j^(E„n - LR„ - MR^) + 2(^—j^{Epn - AfR„ - iVR^) 

+2(t(;iQ,Raa + W2/3R/3/3) - {Wia + ^^2/?) (2-ffRa X R^) 

+2— (R„a - R/3/3) X R,;3 + — " (AC/" + "^(^^ + + N^)]ry 

+5(Vfe, V%, V'is, V'ijj, 4, «|, VJi VJl, 012) . 
On the other hand, we can get by (|5.5ip - (|5.53|) that 

E^a + Ej ME^ MEp {L-N)Ep {L-N)E^ 

- " ~ + 2E 2E ^° 

+5(4,4,V(5^,v4,ai2), 



and by (^5l3\\ and (fOTT) . 



2('(i;iQ,Ra„ + W2l3'Rf3l3) " ("i^la + 'W2l3){2HKa X R/3) 

U'"'(L — N) 
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Thus, we arrive at 



(a-i)(r-r; 



Al^—jEn - 2C/'"(F„R„ + H^K^) + AwiK^ + Au-sR/? 
+2i^—)jEan - LKa - MRp) + 2(— j^(^^n - MR„ - NK^) 

" -^^^ ~ — + — E — E — 

(^R. - ^R. + (L - N)n) + ^ (^R. + ^R, + Mn) 
A?7" + — + 2M2 + iV2)jn + 2^//(^— j^R„ + 2^/?(^— j^R^j 
+5(<5r, (^jf , V(^ij, V6j^, V^6r, V^6^, 6]^, 4, V^i V4, ais) 

r TTi _l_ 

= ^ f 2(M2 - LiV) - Ai? + ^i— 
E \ E 

+ ( - 2U^H^ + Au;i - (L - iV)(— - 2M(— ) Jr„ 

+ ( - 2C/"i?;3 + Azz;2 - 2M(— )^ + (L - iV)(_)^)R^ 
+5(<5r, (^jf , V6r, V6j^, V^6r, V^6^, 4, J|, V^l^, V^, ais) . 
And thanks to ()5.3ip , we have 

Aw2 - ' \ / V y 



which together with Lemma 15.61 imphes that 
(A-1)(R-R)i 

= ^{2iM^-LN)-AE+^^^)n 

Tjn Tin 

+—{L + N- 2EH)ji^ + —{L + N- 2EH)^R0 

2 

+ ^(V'^'^fl, V^5^, V'^.^i V'^aia) 

fc=0 

2 

= Yl V'=<^;j, V'^Ji V'^ais) . (5.56) 

k=0 

Now we are position to prove (|5.32|) . Firstly, by (|5.2p and the fact that the initial surface 
is parameterized by the isothermal coordinates and (j3.1ip holds for t = 0, we know that all 
the relations in (|5.32p hold for t = 0. Hence, 

5a(0) = 0, 5i2(0) = 0, 5^(0) = 0. 
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Taking the energy estimate to (|5.33|) and (|5.34p . we obtain 

||5a(t)||L2 <C[ (||5a||L2 + l|ai2||L2 + ¥R\\m + ¥R\\m)dT, 

Jo 

l|ai2(i)llL2 <C {\\6a\\L2 + ||ai2||L2 + \\6ji\\h2 + \\6j^\\H2)dT. 

Jo 

Using the elhptic estimate, we deduce from (|5.35|) - (|5.38p that 

\\6n\\H2<C\\6j^\\Hi, 

012)11^2 < C(||5i2||l2 + \\5r\\h^ + \\S^\\hi), 

and from (j5.56p . it follows that 

\\Sjlit)\\H2 <C I (||(5a||L2 + ||Si2||l2 + WSrWh-^ + ||5jj||H2 + || (5^;, <5|, 012) ||//2) dr. 

Jo 

Thus, we obtain 

¥a{t)\\L2 + \\ai2{t)\\L^ + \\5^{t)\\H2 

<C [ {\\6a{T)\y + ||5i2(t)|U2 + \\6j^{T)\\H2)dT, 

Jo 

which implies ()5.32p by Gronwall's inequality. 
5.5 Remark on the general case 

In this subsection, we describe how to adapt our method to deal with the case in which the 
surface is parameterized by a finite number of isothermal coordinates. Assume that we need 
A'" local chart to parameterize the initial surface So = U^L^Sq where each Sq is open and 
parameterized by isothermal coordinates: 

Rj,(xi,X2) : Jl' ^ 5^, l<i<N. 

Let {ip^}i<i<N be a partition of the unit subordinate to {5g}i<i<Ar; that is, 

N 
i=l 

At each local chart, R* is defined by 



dt 

where W2) is defined by 



'E- R\y^,L E, ■ 



33 



While, {v, H, 11) is determined by solving the following system: 
v% = 2Hv'', 

2— = a'^'v'^^^ + v^b'lhi + 27;"F„, 

see Section 2 for some notations. As the above equations are coordinate- invariant, (f , 11) 
does not depend on the choice of coordinates. In this case, the energy functional is given by 

£{t) = ||(-A5jV||i.(5^) + \\{-As.fvm,^s,) + \\{-^s.?H\\l,^s,y 

where is the tangential component of the velocity, and is the Laplace-Beltrami op- 
erator on the surface St at time t. In the isothermal coordinates, A^j = -^A. Then, as in 
section 4, we can obtain a uniform estimate for £{t). Let {4''^{t^xi^X2)}i<i<N be a partition 
of the unit on St given by 

</>^(t,xi,X2) = V'H;R^(^i^^2)) ;^\t,X)=i:\R\,o{R\t)r\X)). 
Ej=i V'n^>^H*>3;i,2;2)) 

Indeed, we have 

N 



i=l * 

^ r 1 1 

y / (P\-—A)''v^i-—A)''+^HEidxidx2 + L.W.T. 
-, J Ei bji 

[ {-As,)''v''{-As,)''+'HdSt + L.W.T.. 

JSt 



And, similarly. 
d 



V / (t)\-—^)''H{-—Af+^v''Eidxidx2 + L.W.T. 
I {-As,fH{-Asy+^v''dSt + L.W.T., 

JSt 



where L.W.T. denotes the lower-order terms. Thus, we have 

|(||(-A5jV||i.(<,^) + \\{-As,)^H\\l,^St)) = L-W.T. 
6 Appendix 

6.1 Derivations of the equation ( 12.141) and the energy law 

In this subsection, we give the derivations of the equation (j2.14p and the energy law in the 
case Bai3 = Baap. The reader can also find a short version of the derivation in |10j . 
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Firstly, we have 

^ cf^^iBa^s - b^s) = 2{kiB - {h - ei)H)a^(^ + 2eib"^. 
We get by (H^D that 

= M-%^^^p + {M-%b^p + Ml^^)^. 
Using = 6a7,/3, we infer that 

= = ^^''K^c. = ^^'K,, = 2a-^H^, 
^a-^{b^'b^s),a = a^^b^'b^s,a = a^^b^%^,s = b^'b^^, = 6"^6^,„, 
from which, we can deduce that 

M^%^ = 2{hB-{h-ei)H)b'^g + 2eib^^'bl^ 

= -AkiHa'^^B^a - a"^(2(A;i - ei)H^ + eib^^b^s - ^kiHB) ^, 
M^f^ = 2kia''^B^^p-2{ki+ei)a^^H^^p. 
Let K be the Gaussian curvature. It is easy to see that 

b^'^bo.p = - 2K, b'''^b]^bay = 2H{AH^ - 3K), 

which imphes that 

M'^'^b^^bap = Aki{B - H){2H'^ - K) - 8eiH{H'^ - K). 

Let P =^ n + 2{ki - ei)H'^ + 2ei{2H^ - K) - AkiHB. Then we obtain 

= -{Ua'^^a^),^ + 2eo(5"^a^),„ + M'^^^s.p + (M'^'^bf^bo^p + Mf^)n 
= (Pa"^a^),« + 2eo(S"^a;3),„ - Ak^H a^^ B ^^Sip 

+ (^2A;i(ArS - 2KB) - 2{ki + ei)(Ari? + 2H{H'^ - K))^n. 

We stih use 11 to denote P. Thus, (|2.14p fohows easily. 

Now we derive the energy law of ()2.14p . We infer from ()2.14p that 

= / V • ( - (na^/^a,),^ + 2eo(5"^a,),;3 - AkiHa'^^ B^p^^ 

+2ki {ArB + B{b'^b^^ - AH"^)) n - 2/ii {ArH + H{b'^b^ - 2H^)) nj dS 

= y na° • - 2eoS"^a„ • v,^ - AkiHa'^^ B^/^Va 

+2ki{AYB + B{bpi - 4^2)) - 2fii{ArH + H{b'^b^ - 2H^))v'^dS 
= j -2eoS'^PSo,p + 2kiB(2{Hv'')^^ + {bpi - AH^)v'' + Apz;") 

-2iiiH[Atv'^ + {bpi - 2F2)t;")d5, ^1 = A;i + ei. 
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Here we used — 2Hv^ = such that 



On the other hand, the Helfrich energy can be simphfied as 

Eh = j ikiiH - B)^ + 4:ei{H^ - K)dS, 

where K = ^(4i?^ — b'^ba) is the Gaussian curvature. As J KdS is a constant independent 
of the time, we have 

-Eh = J^{8hiH-B) + 8eiH)—dS 

= aJ ifiiH - kiB){Arv'' + v'^bpi + 2v''H^^). (6.2) 

Adding up ()6.ip and ()6.2p . and using v'^^ — 2Hv^ = again, we obtain the fohowing energy 
dissipation law: 



1 d 
2dt 



+ j ^>|vpd5) = -2eo j S^^SapdS. 



6.2 Some basic estimates in Sobolev spaces 

Let us first recah some product estimates and commutator estimates. 

Lemma 6.1 Let s > 0. Then for any multi-index a,/3, it holds that 

< C(||/||L°°||fl'||iys+|c| + |/3| + ||fl'||L°°||/||H''+l<»l + l^l)- 

In particular, we have 

\\fg\\Hs<C{\\f\\L^\\g\\H^ + \\g\\L^\\f\\Hs). 
Lemma 6.2 Let s > and F(-) G C°°(R+) with F{0) = 0. Then 

\\F{f)\\H^<CmL-)\\f\\H^- 

Lemma 6.3 Let s > 0. It holds that 

\\[A',9]f\\L^<C{\\Vg\\L^\\f\\Hs-i + \\g\\Hs\\f\\L^). 

Here A = (-A)3. 

Lemmas 16. 1116. 31 are well-known, see |1H [M] for example. 
Lemma 6.4 Let s > and k > 1 be an integer. Then it holds that 

||(aA)'=/||^,. < C{\\a\\H2){\\f\\Hs+2, + \\a\\Hs+24f\\m). 
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Proof. We will prove it by induction on k. For A; = 1, using Lemma l6.ll and the Sobolev 
inequality, we get 

||aA/||j^. < C(||a||i^||/||^.+2 + ||a||^,s+2||/|Uoo) 
< C{\\a\\H2\\f\\Hs^2 + \\a\\Hs+4f\\m)- 
Assume Lemma 16.41 holds for /c — 1. Then using the induction assumption, we have 

||(aA)V||L2 = \\{a^f-\aM)\\L^ 

< C(||a||j:f2) (||aA/||^s+2fe- 2 + II a|| j:^s+2fc-2 ||aA/ 11^:^2) . 

We get by Lemma |6. II that 

\\aAf\\Hs+2k-2 < C(||a||^2||/||j:^s+2fe + ||a||j^s+2fe||/||j^2), 

||a||j|^s+2fc-2||aA/||^2 < C||a||j|^.+2fe-2(||a||^4||/||^2 + ||a||^i-2 1|/||^4) 
< C'(||a||^2)(||/||^.+2. + ||a||^.+2.||/||^,2). 

Here we used the following interpolation inequality in the last inequality: 

||«||h4 < ||o||//2 ||a||j:^s+2fc ) ll'3^llH''+2fe-2 ^ ll'3^ll//2 1 1 « 1 1 j:/3+2fe , 

with e = {s + 2k- 4)/(s + 2k- 2). Thus, we get 

||(aA)'=/||H= < C(||a||^2)(||/||^.+2. + ||a||^.+2.||/||H2). 

The proof is completed. □ 

Lemma 6.5 Let s > and sq G (1i2). Assume that a > cq for some positive constant cq. 
Then we have 

||(aA)V||H» >c||/||h.+2. -C(||a||^.o+i)ll«llH»+2.||/||^^.o. 
Proof. We prove the lemma based on the induction assumption on k. For fc = 1, we have 
||(aA)V||H» >co|||A^+'||L2-||[A^a]A/||i2. 

We write 

[A^ a] Af = [A^A, a]f- 2A^(Va • V/) - A^(Aa/), 

which along with Lemma 16.11 and Lemma 16.31 implies that 

||[A^a]A/||i2 < C(||a||^.o+i||/||^.+i + ||a||^,.+2 ||/||^^.o). 

This yields the case of A; = 1 by an interpolation argument. 

Now let us assume that Lemma 16.51 holds for k — 1. Using the induction assumption, we 
have 

||(aA)VllL^ = ||(aA)'=-i(aA/)||i2 

> c||aA/||j|^.+2fc-2 - C(||a||j|^so+i)||a||j;^.+2(fe-i) ||aA/||/i-.. 
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Using the case of A; = 1, we get 

ll<^^/llj?'>+2'=-2 > C\\f\\fj2k+s - C{\\a\\Hsa + l)\\a\\fjs+2k\\J \\H=0 , 

and by Lemma |6.H 

||a||j:^s+2(fe-l) ||oA/||j7s 

< C'||a||j;^s+2(fe-l) (||a||j:^so + 2||/||/fso + 1 1 « 1 1 //''O 1 1 / 1 1 //^o + Z ) 

< e||/||^.+2fc +C(||a||^.o+i)l|a|l/f'>+2'=ll/llH-o. 

Here we used the following interpolation inequality in the last inequality: 

.... II 110 Mill 6 MM Mill — 011 110 

^ Ipll/f^O j:^s+2fc 5 |p||/i'«+2'=-2 < 1 1 ^ 1 1 /f^o 1 1 ^ 1 1 j|^s+2fe , 



||(aA)Vlk^ >c||/|| C{\\a\\ 



with = 2/{s + 2k — sq). Taking e to be small enough, we obtain 

The proof is completed. □ 
Lemma 6.6 Let s > and k > 1 be an integer. Then we have 



[ai,(aA)'=]/||^,. <C(||(G,ata)||^2)(||/||^.+2. + ||(a,5ta)||^.+2.||/||/^2^ 
[V,iaA)'^]f\\Hs < C7(||a||^3)(||/||^.+2. + ||a||H.+2.+i ||/||^,2) . 



Proof. We write 

k-l 



[dt,{aA)'^]f = ^(aA)^[9i,aA](aA)'=-^-V 
= ^(aA)^(5iaA)(GA)'=-^-V. 



Then the first inequality can be seen easily from the proof of Lemma 16.41 The proof of the 
second inequality is similar. □ 

Lemma 6.7 Let s > and k > 1 be an integer. Then we have 

Proof. As in Lemma [6. 5 ^ this lemma can be proved by the induction argument, however, we 
omit the details here. □ 



38 



6.3 Elliptic estimates 

We consider the following elliptic system: 



(6.3) 



Wi W2 



We write 

I 

Then (16. 3p is reduced to solve the following Poisson equations: 

A0 = /i, Ai, = h. 

Thus, we have 

Lemma 6.8 Let s > I. If E e H%T'^),fi,f2 G H'-^T"^), then the system (E^) has a 
solution (Wi,W2) satisfying 

\\Wi\\hs-i + \\W2\\hs-i < C{\\E\\Hs){\\fl\\Hs-i + ||/2||//-i)- 

Next we consider the elliptic equation: 

- AU + aU = f. (6.4) 

Lemma 6.9 Let s > 0. Assume that f G H''{T'^), and a € n L°°{T'^) with 

a{x) > 0, / a{x)dx > oq > 0. (6.5) 

Then there exists a unique solution U € ff*"''^(T^) to |g.^[ ) satisfying 

\\U\\Hs+2<C\\f\\H^. 
Here C is a constant depending only on and ||a||_f/TiL°°- 

Proof. The proof of existence part is standard. Here we only prove the estimate. Taking the 
L^ inner estimate gives 



IVC/III2+ [ a\U\'^dx= [ fUdx. 



Let U = /rp2 Udx. Therefore, we have 

IIC/III2 < 2\\U - tJ\\l2 + 2||i7||^2 < 2||VC/||i2 + 87r2|i7|. 
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On the other hand, we have by (j6.5p that 



ao|C/p < / a\U\'^dx < 2/ a\U\'^dx+[ a\U-U\'^dx 

< 2 / a|;7p(ix + ||a||Loo||VC/|||2. 

This yields that 

\\U\\l2 <C{\\VU\\l2+ [ a\U\'^dx). (6.6) 

Using the elliptic estimate in H^, we obtain 

\\U\\hs+2 < C{\\aU\\H^ + ll/ll/f.) < C(||a||Loc,||C/||^,. + \\a\\Hs\\U\\L^ + \\f\\Hs), 

from which and (j6.6p . the desired estimate follows from an interpolation argument. □ 
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